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Abstract 


This paper introduces a method for the evaluation of the occurrence of tipping 
points based on the combination of probability intervals from experts judgments 
elicited in face-to-face interviews. The computation of such conditional probabil- 
ities is based on the aggregation of imprecise probability judgments through the 
Steiner point. The probability of a tipping point can be updated by the standard 
Bayes rule to generate tipping point scenarios. Our results suggest that tipping 
events may happen with relatively large probabilities, in contrast with the view 
that tipping points are low-probability-high-impact events. This suggests that 
mitigation and containment policies cannot be further postponed. 


Keywords: Bayesian updating; aggregation of opinions; global warming; judgmental 
forecasting; Steiner point; tipping points 
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1 Introduction 


As of today, there is a large literature suggesting that the consequences of climate 
change may involve abrupt changes and tipping points (e.g. Lenton et al., 2019).! 
Armstrong McKay et al. (2022) identify a finite core of tipping point elements able to 
modify the Earth system functioning and show that six climate tipping points global 
elements are likely even in the Paris Agreement range of 1.5°C to 2°C warming. The 
concept of tipping point was introduced into the scientific debate in the '80s of the last 
century by the IPCC to represent large scale discontinuities in the climate system. At 
that time, experts believed that tipping points would be crossed if global warming had 
exceeded 5*C. In the 2018 IPCC report (IPCC, 2018b) experts suggested that tipping 
points could be crossed even between 1 and 2*C of global warming.? 

The main objective of this paper is to introduce a new methodology for the 
evaluation of occurrence of tipping points based on the combination of probability 
intervals from experts judgments elicited in face-to-face interviews. When experts face 
uncertainty or have imprecise knowledge about future states of the world, alterna- 
tive approaches that differ from Bayesian pool methods are required. In this paper 
we apply an aggregation method based on the Steiner point introduced in Basili and 
Chateauneuf (2020). This method assumes that experts have imprecise information 
represented by convex sets of probability distributions with the requirement that the 
intersection of these sets is not empty and experts have at least one common probabil- 
ity distribution. A desirable feature of this approach is that the opinion pooling is the 
probability distribution of the tipping points we considered. This differs from other 
methods where some type of (simple or weighted) average is defined (see, e.g., Clemen 


& Winkler, 1999, 2007). Furthermore, by means of a standard Bayesian updating rule 


1A tipping point is defined as “[a] level of change in system properties beyond which a system reorganizes, 
often abruptly, and does not return to the initial state even if the drivers of the change are abated. For 
the climate system, it refers to a critical threshold when global or regional climate changes from one stable 
state to another stable state” (IPCC, 2018b). See also the recent article by Armstrong McKay et al. (2022) 
for a review and a thorough definition of tipping point. 

2Lenton et al. mention that some “models suggest that the Greenland Ice sheet could be doomed at 1.5°C 
of warming, which could happen as soon as 2030” (Lenton et al., 2019, p. 592). 
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we can provide posterior probabilities that may be used to generate future scenarios 
or hypotheses for future studies, once more data become available (see Section 3 for 
more details on the aggregation approach).? It is important to notice that the con- 
struction of the Steiner point requires that the intersection of the experts’ opinion be 
non empty. However, if that were not the case, the aggregation could be applied to 
subsets of consistent experts (see Section 3.1). 

To provide an assessment of the occurrence of tipping points we use the elicitation 
data in Kriegler et al. (2009).* In Kriegler et al.'s seminal paper, experts are asked 
to provide probability intervals about three temperature scenarios. In this context, 
intervals of probabilities are a representation of uncertain and imprecise judgements.? 
The aggregation method proposed in this paper is new in the context of problems 
related to climate change. 

The paper proceeds as follows. Section 2 briefly discusses the related literature. 
Section 3 introduces the fundamental aspects of our theoretical framework. In Section 
4 we present the conditional probabilities of occurrence of various tipping points com- 
puted via the Steiner point as well as the results of the Bayesian updating. Finally, 


Section 5 concludes the study. 


2 Related Literature 


Communicating uncertainty about climate change is crucial to effectively influence 
policy decisions and shape public opinion. In this context, the IPCC special report 


Global Warming of 1.5°C (IPCC, 20182) is an important example. The IPCC special 


3It seems that there are no follow up studies on experts” assessment of tipping points. This is also 
confirmed in the correspondence with leading researchers in the field. Other recent studies use the data and 
results in Kriegler, Hall, Held, Dawson, and Schellnhuber (2009), see, e.g., Gaucherel and Moron (2017); 
Wunderling, Donges, Kurths, and Winkelmann (2021). 

*The data are collected from Kriegler et al. (2009) supplemental appendix and are available at this link 
along with the replication code. 

5In Kriegler et al. (2009) the low temperature scenario (Low) considers an increase between about 1*C 
and 2*C by year 2200 in comparison with year 2000. The medium temperature scenario (Medium) considers 
an increase between about 2°C and 4°C by year 2200 in comparison with year 2000. The high temperature 
scenario (High) considers an increase between about 4°C and 8*C by year 2200 in comparison with year 
2000. 
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report is based “on the assessment of around 6,000 peer-review publications, most 
of them published in the last few years” (IPCC, 2018a, p. v). The report aggregates 
multiple forms of knowledge to address and communicate the degree of certainty (or 
lack thereof) of specific findings. 

In general, however, some critical issues emerge in the treatment of uncertainty: 
findings are based on multiple lines of evidence and are expressed using confidence 
qualifiers and many of them depend on certain model assumptions. In addition to 
that, findings have to be updated if new information becomes available. Borsuk and 
Tomassini (2005), Tomassini, Reichert, Knutti, Stocker, and Borsuk (2007), Knutti et 
al. (2008), Zickfeld et al. (2007), Zickfeld, Morgan, Frame, and Keith (2010), Kriegler 
et al. (2009) highlight that aggregation of probabilistic projections with a variety of 
statistical models is an unsolved problem. One fundamental challenge of the assessment 
process is to summarize such information into a unique quantity. However, due to the 
heterogeneity in sources and quality of the information, obtaining a unique synthetic 
measure may become a daunting task. 

When formal statistical procedures are unavailable, expert judgment approaches 
are often employed to provide an assessment of uncertainty (see e.g. Mastrandrea et 
al., 2011).* Significantly, “one option is to resort to imprecise probability (e.g., Kriegler 
and Held (2005); Hall, Fu, and Lawry (2007); Tomassini et al. (2007)), that is, consider 
an uncertainty in PDFs [probability distribution functions] or sets of PDFs” (Knutti 
et al., 2008, p. 2658). If each expert has a set of probability distributions, the mean 
value for each scenario is evaluated along with all the individual PDFs (Tomassini et 
al., 2007). 

Experts” quantitative judgments are often elicited in face-to-face interviews. Then, 
using a range of different procedures, mean (averaged over experts) ranks are computed 


(Kriegler et al., 2009; Zickfeld et al., 2010). In Zickfeld et al. (2010), 14 experts (leading 


SIn 1975, the U.S. Nuclear Regulatory Commission (NCR) introduced for the first time a procedure for the 
elicitation process and, since then, techniques and methods have spread to other areas such as volcanology, 
public health, ecology, aeronautics, climatology etc. (Cooke, 2013). 
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climate scientists) discuss about three scenarios (high, medium, low) of net radiative 
forcing at the top of the atmosphere from anthropogenic sources through the year 
2200. Experts use a cardinal scale from 0 (no chance) to 1 (definite chance) for each 
of the three forcing trajectories. Experts elicit probabilities for each scenario and in 
the following step they are asked to estimate the median of the mean trajectories of 
warming between 2000 and 2050. What they find is that it falls between 0.16*C/decade 
and 0.36°C /decade. 

Providing a realistic assessment of the probability of a tipping event may have a 
fundamental impact in influencing policy decisions concerning climate change. Unfor- 
tunately, there are very few models able to simulate abrupt climate changes as 
consequences of realistic external forcing and, even when sophisticated climate mod- 
els show effects, in particular a temperature response, they are weak with respect to 
Dansgaard-Oeschger events.” To manage the complexity and the uncertainty implied 
in climate models, policy makers (PMs) have no alternative but to resort to experts 
and obtain experts’ elicitations; that is, expert probabilistic judgments (Colson & 
Cooke, 2018), that are represented by imprecise probabilities or intervals of probabil- 
ities, generally. Once the elicitation is completed, the data are aggregated. Using a 
reliable aggregation method to provide a synthesis of the various opinions seems of 


paramount importance, also given the fact that no clear dominant approach exists.? 


"Dansgaard-Oeschger events are rapid climate fluctuations such as the Greenland ice melting occurred in 
the Eemian interglacial. This event took the form of rapid warming episodes, followed by gradual cooling 
periods, that increased the average annual temperature on the Greenland ice sheet of 8 °C over 40 years 
(Dansgaard et al., 1993; Heinrich, 1988). 

$See Lam and Majszak (2022) for a recent review on the challenges and opportunities of expert judgment 
for the assessment of climate tipping points. 
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3 Theoretical Framework 


The aggregation process of probabilistic opinions (opinion pooling) entails a function, 
known as a pooling rule, that elicits a consensus distribution. Such a consensus distri- 
bution is determined among not necessarily independent and fully competent experts 
when each of them has multiple priors on future states of the world. 

Under uncertainty or deep uncertainty experts have partial, incomplete or fuzzy 
knowledge and their beliefs cannot be represented by a unique, additive and fully 
reliable probability distribution, but either by a finite set of them, an interval of 
probabilities or by a non necessarily additive measure (e.g. a capacity) (e.g. Basili & 
Chateauneuf, 2011, 2016; Basili & Pratelli, 2015, and references therein).? 

In this paper we assume that opinions are expressed through different probabil- 
ity distributions and that there exists a PM who adopts a multiple priors decision 
model. The set of probability distributions of all experts can be considered a reflection 
of the PM's assessment of the reliability of available information about the underly- 
ing uncertainty, that is, her perception of uncertainty; the elicited aggregation rule 
incorporates the PM's attitude about scanty and vague information. Facing the set 
of all probability distributions attached by experts to possible events, the PM evalu- 
ates their probability intersection. Such a common opinion is the Steiner point of the 
convex capacity that emerges from the aggregation of experts' opinions. 

The aggregation via the Steiner point provides a natural and tractable way to 
summarize opinions and only assumes that experts have imprecise information but 
nothing about their competence, experience and independence. This is, for example, 


the case of non Bayesian groups, i.e. groups where experts do not have a unique and 


“Deep uncertainty includes “situations in which we are still able (or assume) to bound the future around 
many possible plausible futures and situations in which we only know that we do not know” (Marchau, 
Walker, Bloemen, & Popper, 2019). See also Rohmer, Le Cozannet, and Manceau (2019) and Frederikse et 
al. (2020) for some recent discussion of the concepts of uncertainty in the context of climate change issues. 
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additive probability distribution.!? In addition, the Steiner point not only allows us to 
derive the consensus distribution about a given event, but it is amenable to Bayesian 
updating. This means that the PM can simply update the elicited distribution when 
new information is made available without calling for a new round of interviews. 
Whenever experts have opinions that are too heterogeneous the aggregation cannot 
be performed (see the G-consistency condition below in this section). Yet, the Steiner 
point can be applied to subsets of individuals allowing us to identify different trends 
and common opinions among consistent experts. In general, the aggregation through 
the Steiner point suggests that the experts are consistent and that the interpretation 
of a tipping point as a high-impact-low-probability event seems to be incorrect (see 
also Lenton & Ciscar, 2013; Lenton et al., 2019). We find, in fact, that the probability 
of crossing a tipping point threshold is non trivially larger than zero in nearly all cases, 
even for low climate change scenarios. Furthermore, it has been claimed (Lenton & 
Ciscar, 2013) that a more adequate picture for the representation of tipping points 


would be to provide the (joint) probability distribution for tipping each element. 


3.1 Experts’ Consistency 


Let us consider a finite set S = (51,...,5,] of states of the world and let X = 23 be 
the o-algebra associated to the set S. Since the study deals with discrete events, P is 
a probability mass function (pmf) on (S,2) such that P : X > [0,1]. For a specific 
pmf, p(s;) = P(S = si) for any i € {1,...,n} and 35; 4 p(si) = 1. 

Consider now a non-negative set function v : 3 — R. Such a function is a capacity 
(a non additive probability) on (S,%) if, for any A,B € Y, ACB = v(A) < v(B), 


v(O) = 0 and v(S) = 1, where Ø is the empty set. The capacity v(-) is convex if 


10 The experts that provide judgements for the IPCC assessment may be interpreted as a non Bayesian 
group, as they typically provide evaluations in terms of qualitative levels of confidence and quantitative 
likelihoods of occurrence of a given event generally expressed as intervals of probabilities. Mastrandrea 
et al. (2011) characterize the degree of certainty in key findings with two metrics: "[c]onfidence in the 
validity of a finding, based on the type, amount, quality, and consistency of evidence (e.g., mechanistic 
understanding, theory, data, models, expert judgment) and the degree of agreement, [c]onfidence is expressed 
qualitatively" and "[q]uantified measures of uncertainty in a finding expressed probabilistically (based on 
statistical analysis of observations or model results, or expert judgment)". 


we V(AUB)+v(ANB) > v(A) + v(B). The corresponding dual capacity is defined as 
159 (A) = 1 — v(.A^). The dual capacity v(.A) is concave. Uncertainty is modeled via the 
: core of the convex capacity v. The core C(v) is a set of probability distributions P on 
11 . ($,X) such that P(A) > v(A) VA € X. 

162 The PM's problem is to learn the distribution of a certain event by aggregating 
13 experts’ opinions. Let Po denote the unobserved pmf that governs the phenomenon 


1 under study. The PM asks expert j, (j = 1,...,m), to provide a lower and upper 


1s bound for the probability p; = po(s;) = Po(S = s;). The set of possible probabilities 


16 considered by expert j is defined as 


127 The PM will accept the expert’s opinion if two consistency conditions are met. Specif- 


we ically, for the expert's opinions to make sense the bounds a? and b must meet the 


19 following conditions of individual consistency (I-consistency). 


mo Condition 1. | I-consistency] The bounds a? and bÍ satisfy the conditions 0 < 
m al <b Sand 5) al <1< y, b. 


172 The set PÍ is not empty if and only if the I-consistency condition is met. 


13 Furthermore, PJ can be seen as the core C(v/) of a given convex capacity v/ where 


v (A) = max b» al, 1— >` b (1) 

i€(is;c.A) i€ (isi. A) 
ia and the set A C S is a set of states of the world (Chateauneuf & Cornet, 2018; 
vs De Campos, Huete, & Moral, 1994). It is expected that the unknown distribution Po 


ve be in the the set Pi of the generic expert j, ie. Py € PÍ. In addition to that, it 
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is expected also that P) be in the intersection of all P7. Hence, the following group 


consistency (G-consistency) condition is supposed to be met. 


Condition 2. [ G-consistency] The intersection of the probability sets associated 


with the pool of experts is non empty, i.e., P = ne, Ps = © and that Po € P. 


Whenever the intersection set 7 = O, this is, when experts have conflicting 
opinions, the PM may still be able to extract valuable information by applying the G- 
consistency condition to a subset of experts. This approach may reveal, for example, 


whether some experts are more or less optimistic about a given phenomenon. 


3.2 Aggregation and Scenarios 


The consensus distribution set is associated to a convex capacity v(-), defined as 


v(A) = max 5 ai, 1— X bi (2) 


1€[1:s¡E A} 1€Efi:s¡ ZA} 


where a; = max; al and b; = min; bi . In this context the Steiner point is particularly 
relevant, since the Steiner point of C(v) is the center of the core of the convex capacity 
v and represents the consensus probability for the given set of experts. The Steiner 
point is denoted as II" € C(v). 

It is interesting to notice that when the set of states of the world is finite, the 
Steiner point coincides with the Shapley value (Basili & Chateauneuf, 2020; Pechersky, 


2015; Shapley, 1971) and it is easily computed via the following expression 


r= x Meter la HA — o(A\{s;}), i= 1,..., n. (3) 


l 
n: 
s¡EACS 


194 


195 


196 


197 


198 


199 


200 


201 


202 


203 


204 


205 


206 


207 


208 


Equation (3) shows how the Shapley value represents the average marginal individual 
contribution over all the possible different permutations in which the grand coalition 
S may be formed (Basili & Chateauneuf, 2020).!! 

Our problem consists of two states of the world (n = 2), this is, whether a tipping 
point occurs (Tip) or it does not occur (No Tip). The two states of the world occur 
with, say, probability II? and I = 1— IIT, respectively. Let us now define the random 
variables T € [Tip, No Tip) denoting the occurrence (or not) of a tipping point and 
C € (Low, Medium, High) representing possible temperature scenarios as defined in 


Section 1. Hence, we have the following Bayes rule 


P(C|T — Tip)P(T — Tip) 
(C|T = No Tip)P(T = No Tip) + P(C|T = Tip)P(T = Tip) 


(4) 


P(T = TipO) = 5 


The Bayes rule in equation 4 can be used to update the probability of occurrence of a 
tipping point obtained via the Steiner point. Specifically, we choose II? for P(T = Tip) 
as a prior probability, while P(C|T = Tip) and P(C|T = No Tip) can be obtained 
from expert knowledge. Alternatively, we can feed the formula a grid of values. In this 
last case the posterior probability P(T' = Tip|C) can be graphically represented as a 


surface. For the empirical analysis in Section 4 we opt for the latter. 


“Recent research suggests that, for independent inputs, the Steiner point-Shapley value is bracketed 
between two different Sobol’ indices. This result seems to hold also for the case of dependent inputs or 
expert judgments (Song, Nelson, & Staum, 2016, see also Owen and Prieur (2017) for further uses of the 
Shapley value in the context of ANOVA). Sobol’ index provides a measure of the importance of inputs to 
a function and is defined in terms of the functional analysis variance decomposition (Sobol’, 1990, 1993). 
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4 Eliciting Probabilities of Tipping Points by the 


Steiner Point 


Kriegler et al. (2009) elicited beliefs of experts about the probability of triggering 
major changes in the Earth system associated to seven tipping points (Table 1 reports 
the tipping points under analysis) for three different global median temperature sce- 
narios. As previously mentioned in Section 1, the three scenarios consider the crossing 
of a tipping point for a given temperature scenario and a fixed time horizon. Specifi- 
cally, Low refers to an increase between about 1°C and 2°C by year 2200 in comparison 
with year 2000, Medium refers to an increase between about 2°C and 4°C by year 2200 
in comparison with year 2000 and High refers to an increase between about 4°C and 
8°C by year 2200 in comparison with year 2000. In light of the new evidence produced 
in over a decade, Armstrong McKay et al. (2022) classify new and known tipping ele- 
ments as tipping points that act on a global scale or only locally (such as the dieback 
of Boreal forests, BOFO in Table 1). For other phenomena such as El Niño (NINO in 
Table 1), there seems to be insufficient evidence for it to be characterized by a tipping 
element. The decline of the ocean carbon sink (DOCS in Table 1), which was previ- 
ously considered a tipping element, is categorized as threshold-free. For consistency of 
exposition with respect to Kriegler et al.’s data, we consider NINO and DOCS as if 
they were tipping elements.!? 

The raw data collected by Kriegler et al. (2009) are summarized in the dumbbell 
plots in Figures A1 to A7 in Appendix A. Each figure consists of three plots and each 
plot refers to one of the three temperature scenarios. To every row in the plots we 
attach one single expert. By direct inspection, we note, as also stressed in Kriegler et al. 
(2009), the tendency of the experts to place high probability in the high temperature 


scenarios. '? 


12Kriegler et al. (2009) report that BOFO and DOCS “were judged by experts to be of more speculative 
nature” 

13 The data are collected from Kriegler et al. (2009) supplemental appendix and are available at Federico 
Crudu's personal webpage (link) along with the replication code. 
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AMAZ | Dieback of the Amazon rainforest. 

BOFO Dieback of Boreal forests. 

AMOC | Reorganization of the Atlantic meridional overturning circulation. 
DAIS Disintegration of the West Antarctic ice sheet. 

MGIS Melt of the Greenland ice sheet. 

DOCS Decline in ocean carbon sink. 

NINO Shift to a more persistent El Nino regime. 


Table 1 Tipping points considered in Kriegler et al. (2009). 


In order to provide policy relevant information Kriegler et al. (2009) combine, 
using an aggregation rule, the probability intervals supplied by the experts for every 
tipping point and every temperature scenario.!^ Their aggregation rule produces upper 
and lower probabilities. Our approach is substantially different. In fact, the Shapley 
value described in equation 3 returns the probability of occurrence of a tipping point 
conditional on a given temperature scenario. 

Figure 1 features the probabilities associated to every tipping point and each 
climate scenario resulting from the aggregation of the experts' opinions via the appli- 
cation of equation 3. In Figure 1 we see that, in general, a high climate change scenario 
(red triangle) is more likely to trigger a tipping point. While this result is not too sur- 
prising, we also observe that there is a certain degree of heterogeneity across scenarios 
and tipping points. Specifically, AMOC has a very low chance to occur for the low cli- 
mate change scenario. All the other tipping points have, for the same scenario a non 
negligible (say, larger than 1096) probability to occur. As we consider the higher cli- 
mate change scenarios, Medium and High, the probability of having a tipping point 
increases. It is remarkable that for High, the probabilities of occurrence for BOFO, 
DAIS, DOCS and MGIS are very high. 

To provide further intuition for the interpretation of the plots in Figure 1, we 
adapt the terminology used in Kriegler et al. (2009). This is, we label as remote the 


prospect of having a tipping point if II] < 0.1, significant if II] > 0.1 and large if 


14Kriegler et al. (2009) introduce a novel aggregation rule called forced consensus pooling. See also Clemen 
and Winkler (1999) and Nau (2002). 
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II; > 0.5. Under these criteria only AMOC for the low temperature scenario has a 
remote chance of being triggered. DAIS and MGIS have a large probability of being 
triggered under the high temperature scenario, while BOFO and DOCS have a large 
probability of being triggered under both the medium and high temperature scenario. 
The remaining scenarios have a significant probability of triggering a tipping point. 
By applying the Bayes rule in equation 4 we update the results obtained via 
the Steiner point to generate tipping point scenarios for different combinations of 
P(C|T = No Tip) and P(C|T — Tip). Every point in the resulting surfaces (figures 2 
to 8) is the conditional probability of crossing a tipping point given the likelihood of the 
(conditional) occurrence of the temperature scenario. We notice that the high temper- 
ature scenario produces large probabilities for the occurrence of tipping points. This 
result is particularly clear for BOFO (figure 3 (c)), AMOC (figure 4 (c)), DAIS (figure 
5 (c)), DOCS (figure 6 (c)) and MGIS (figure 7 (c)). Committing to low climate change 
scenarios may produce remote probabilities of tipping, at least for AMOC (figure 4 
(a)) and DAIS (figure 5 (a)), yet in most of the remaining cases the probabilities of 


triggering a tipping point are generally at least significant. 


13 


1.00- MAZ 1.00- BOFO 
0.75- 0.75- 
E oso E oso 
0.25- 0.25- 
0.00- 0.00- 
0.00 0.25 0.50 0.75 .00 0.00 0.25 0.50 0.75 
V V 
n n 
1.00- AMOC 1.00- DAIS 
0.75- 0.75- 
E" oso E oso 
0.25- 0.25- 
0.00- 0.00- 
0.00 0.25 0.50 0.75 .00 0.00 0.25 0.50 0.75 
V V 
n n 
1.00- DOCS 1.00- MGIS 
0.75- 0.75- 
E oso E oso 
0.25- 0.25- 
0.00- 0.00- 
0.00 0.25 0.50 0.75 .00 0.00 0.25 0.50 0.75 
V V 
n ny 
1.00- NINO 
0.75- 
E“ oso 
0.25- 
0.00- 
0.00 0.25 0.50 0.75 .00 
V 
m 


Fig. 1 Steiner points. The Low, Medium and High climate change scenarios are denoted by a circle, 
square and triangle respectively. II? is the probability that the tipping point is triggered, while I; is 
the probability that the tipping point is not triggered. 
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5 Conclusions 


This paper considers the aggregation of experts’ opinions about uncertain climate 
change scenarios. The aggregation rule based on the Steiner point provides an 
evaluation of the probability of occurrence of a tipping point. 

In this context, the PM can evaluate the individual consistency (I-consistency) 
of an expert, that is, whether an expert is able to provide coherent probabilistic 
evaluations of possible future events, and the consistency of the whole group (G- 
consistency), that is, whether the intersection of the probability sets associated with 
the pool of experts is non empty. If the latter condition is verified the experts in the 
pool share common opinions about future states of the world. The Steiner point allows 
to automatically detect both individual and group consistency, but, more relevantly, 
it can be updated by the Bayes rule, whenever new information is available, without 
repeating a new trial of interviews among experts. 

Ultimately, our results suggests that tipping points have higher probabilities to 
occur under high temperature scenarios. If taken in isolation this result is perhaps 
not surprising. À more interesting and concerning outcome is that the probability of 
crossing a tipping point threshold is noticeably different from zero, even under lower 
climate change scenarios and for the majority of tipping points under investigation. 
The probability surfaces obtained via the Bayes rule show that committing to a low 
climate change scenario may produce remote probabilities of realization of tipping 
points, yet in most cases such probabilities tend to be larger than 1096. Significantly, 
despite the fact that the data are over a decade old, to a large extent the results 
obtained with our methodology overlap with those produced in Armstrong McKay et 
al. (2022), even for the non occurrence of a tipping point for El Niño. 

Due to limited data availability, we can only recover the marginal distribution of 
tipping points. We update the resulting probabilities using a standard Bayes rule for 


all the possible combinations of the likelihood of occurrence of a certain temperature 
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scenario. This approach generates posterior probability surfaces that may be inter- 
preted as tipping points scenarios. The results suggest that some tipping points (the 
dieback of the Boreal forest, the reorganization of the Atlantic meridional overturning 
circulation, the decline of the ocean carbon sink and the melt of the West Antarctic 
and Greenland ice) are extremely likely to occur under the high temperature scenario. 
On the other hand, only committing to low climate change temperature corridors 
would substantially reduce the probability of occurrence of most tipping points. 

Our analysis, consistently with the recent literature on the topic, highlights the 
fact that these type of events may turn out to be one of the main climatic emergencies 
in the upcoming years. Future research can further contribute to the understanding 
of these issues by addressing problems concerning elicitation data, the construction of 
joint distributions for tipping events and the estimation of causal relationships among 


tipping points. 
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Appendix A Figures 


This section displays a set of dumbbell plots that describe the data used in the analysis. 
The plots are similar to those in Figure 1 of Kriegler et al. (2009). Furthermore, in 
Kriegler et al. some of the experts are recognized as core experts. The caption of each 


plot indicates which experts are not core experts. 
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Fig. A1 Elicited probability intervals for the AMOC tipping point. In Kriegler et al. (2009) C9, 
C20 and C22 are not classified as core experts. 
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Fig. A2 Elicited probability intervals for the MGIS tipping point. In Kriegler et al. (2009) M1, M3 
and M14 are not classified as core experts. 
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Fig. A3 Elicited probability intervals for the DAIS tipping point. In Kriegler et al. (2009) D7 and 
D10 are not classified as core experts. 
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Fig. A4 Elicited probability intervals for the AMAZ tipping point. In Kriegler et al. (2009) A1, A6, 
A10 and A12 are not classified as core experts. 
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Fig. A5 Elicited probability intervals for the NINO tipping point. In Kriegler et al. (2009) N3, N6 
and N10 are not classified as core experts. 


B8 B8- B8 

B: B5 B: 
a a a 

£ 

B3 B3- ie 
ú ú ú 

B2 B2- B2 

B1- E B1- - B1- + 

0.00 0.25 0.50 075 00 0.00 025 050 075 00 0.00 025 050 075 
Intervals for Low scenario (BOFO) Intervals for Medium scenario (BOFO) Intervals for High scenario (BOFO) 


(a) Low (b) Medium (c) High 


Fig. A6 Elicited probability intervals for the BOFO tipping point. In Kriegler et al. (2009) B1 is 
not classified as a core experts. 
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Fig. A7 Elicited probability intervals for the DOCS tipping point. In Kriegler et al. (2009) O2 is 
not classified as a core experts. 
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